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ABSTRACT 



Based on drizzled F606W and F814W images, we present quantitative structural 
parameters in the V-band rest-frame for all galaxies with z<l and Isi4(AB)<24.5 
mag in the Hubble Deep Fields North and South. Our structural parameters are 
based on a two-component surface brightness distribution using a Sersic bulge and 
an exponential disc. Detailed simulations and comparisons with previous work are 
presented. The luminosity-size distribution of early-type galaxies is consistent with 
the hypothesis that their structural properties were already in place by z ~1 and have 
evolved passively since then; early-type galaxies were ~ 1.35(±0.1) mag brighter in 
rest-frame V-band luminosity at z^0.7 than now. Compared to present day late- 
type galaxies, those at z~0.7 with Ly>0.2xl0 10 h -2 L show a moderate decrease 
(~30(±10)%) in size (or interpreted differently, a decrease of ~0.77(±0.30) mag in 
the central surface brightness) at a given luminosity. Finally, we make a comparison 
of our results with the infall and hierarchical models. 
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1 INTRODUCTION 

The Hubble Deep Field (HDF) North (Williams et al. 1996) 
and South (Casertano et al. 2000) are among the deepest 
views of the universe ever obtained. These observations were 
specially made to facilitate the study of galaxy evolution. 
The high resolution of the publically available drizzled im- 
ages (0.04"/pix, i.e. ~0.22 h _1 kpc at z^l 1 ) convert them 
into an ideal laboratory for morphological evolutionary stud- 
ies. 

Morphology plays a fundamental role in understanding 
galaxy evolution as each galaxy type could evolve differently 
with z. To facilitate the comparison with the theoretical pre- 
dictions, galaxy evolution studies must be rooted in a quan- 
titative basis. In addition, quantitative morphology enables 
one to study any possible evolution of the galaxy structural 
parameters such as luminosity, size, etc., and their associ- 
ated scaling relations with regards to the values measured 
in nearby samples. 



* E-mail: trujilloOmpia-hd. mpg.de (IT); jalfonso@iac.es (JALA) 
1 We adopt a A-cosmology throughout the paper with Q.M =0.3, 
C A =0.7 and H =100h km s _1 Mpc" 1 . 



In the hierarchical galaxy formation framework, both 
disc and elliptical galaxy properties are expected to evolve 
at high-z. In this paradigm, it is generally assumed that 
much of a galaxy's star formation takes place in a galac- 
tic disc which formed by gas infall into dark matter halos 
(White & Frenk 1991; Kauffmann et al. 1999; Somerville & 
Primack 1999). According to Mo, Mau & White (1998) the 
size of the disc galaxies, at a fixed halo circular velocity, 
must decrease with redshift as R txH~ (z) where H(z) is the 
Hubble parameter at redshift z. On the other hand, E/SO 
galaxies are expected to form via repeated merging (Toomre 
& Toomre 1972). The formation epoch of the E/SO galaxies 
is sensitive to the assumed values of cosmological parame- 
ters, occurring at higher redshift at lower Qm (Kauffmann 
& Chariot 1998). 

There are alternative theories of galaxy formation. One 
is that all galaxies, including E/SO, formed very early via 
monolithic collapse of the gas at high redshift (Eggen, 
Lynden-Bell & Sandage 1962; Larson 1975; Chiosi & Car- 
raro 2002). According to this theory E/SO properties will be 
already in place at high-z. For disc galaxies, some "back- 
ward" approaches have been tried. These approaches use 
the detailed studies of the profiles of the Milky way and 
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other nearby galaxies to propose radially dependent models 
of star formation that could be used to make predictions 
about high-z disc evolution (Cayon, Silk & Chariot 1996; 
Bouwens, Cayon & Silk 1997; Roche et al. 1998). However, 
both the hierarchical and "backward" models predict rela- 
tively similar amounts of evolution in disc global properties 
(size, luminosity, etc.) to 2~1 (Bouwens & Silk 2002). 

The field covered by both HDFs is not large compared 
to others works but they are so deep that they are signif- 
icantly less affected by the incompleteness problem affect- 
ing other less deep samples. Discrepancies associated to in- 
completeness can be found in the literature, for example, 
studying discs galaxies, Schade et al. (1995, 1996) who used 
early ground- and space-based images of galaxies from the 
Canada-France Redshift Survey (CFRS) and Roche et al. 
(1998) in the Medium Deep Survey found a net increase in 
the surface brightness and a net decrease in the size. How- 
ever, Simard et al. (1999), who conducted a detailed analysis 
of the selection effects in the Deep Extragalactic Evolution- 
ary Probe (DEEP) survey found no evolution in disc prop- 
erties up to 

In the following work we conduct a detailed morpholog- 
ical analysis of all the galaxies with z<l in the HDF-North 
and South on the drizzled F606W and F814W images down 
to i8i4(AB)=24.5 mag. Quantitative morphological analysis 
can be conducted using non-model based techniques like the 
concentration and asymmetry of the galactic light (Abraham 
et al. 1996). However, these techniques have the disadvan- 
tage that they do not allow the study of the different struc- 
tural components of the galaxy separately, and particularly, 
the evolutionary analysis of the well known local scale rela- 
tion can not be addressed. 

Improving upon previous analyses, our structural pa- 
rameters are based on a two-component surface brightness 
distribution using a Sersic r 1 ^" bulge (Sersic 1968) and an 
exponential disc. A previous detailed morphological analy- 
sis of the HDF-N was carried by Marleau & Simard (1998, 
hereafter MS98). However, that work only considered a r 1 ' 4 
bulge, instead of a free r 1//n . Our assumption of a n free bulge 
is more in agreement to what it is observed in the nearby 
universe (see, e.g. Graham 2001 and references therein). On 
the other hand, even though the nearly identical quality of 
the HDF-S, a systematic quantitative morphological analy- 
sis of this field have been not yet conducted and it is of great 
interest to compare the galaxy properties between these two 
twin fields. 

In this paper we present the galaxy structural parame- 
ters of both HDF fields and compare our results to previous 
classification schemes. To assess the degree of evolution, if 
any, the early- and late-type luminosity-size relations are 
compared with those observed in the nearby universe pro- 
vided by large local surveys such as the Sloan Digital Sky 
Survey (York et al. 2000). In a second paper (Aguerri & 
Trujillo 2004), we analyse the redshift evolution of the struc- 
tural components of the galaxies (bulge and disc) separately. 
The outline of this paper is as follows: section 2 describes 
the characteristics of the observed fields. Our quantitative 
morphological analysis algorithm is explained in section 3. 
In section 4 we detail our object morphological classifica- 
tion and in section 5 we analyse the accuracy of our analy- 
sis by comparing our results with simulations and previous 
work. Section 6 deals with the issue of selection effects. The 



luminosity-size relations of early- and late-type galaxies are 
presented in section 7. We discuss our results in section 8. 



2 DATA 

Each HDF image consists on a w 2.5 x 2.5 arcmin 2 Wide 
Field Planetary Camera 2 (WFPC2) pointing (Williams 
et al. 1996, HDF-N; Casertano et al. 2000, HDF-S). The 
images were acquired through the broad band F300W, 
F450W, F606W, and F814W filters. For the present work 
we have used the publically available reduced F606W and 
F814W Version 2 Drizzled images with a sampling of 
0.04 arcsec/pixel and a total exposure time of 109050 
s (F606W;HDF-N), 123600 s (F814W;HDF-N), 97200 s 
(F606W;HDF-S) and 112200 s (F814W;HDF-S). The high 
sampling of these images avoids the undersampling problems 
affecting other Hubble Space Telescope (HST) images. 

In order to assure a reliable estimation of the structural 
parameters (see the simulations section) we have selected 
all the l8i4(AB)<24.5 galaxies with z<l. In addition, to 
maintain the analysis in the V-band rest-frame, we have 
also studied the l8i<i(AB)<24.5 galaxies with z<0.4 in the 
F606W filter. The above restrictions leave us with 123 ob- 
jects in the HDF-N and 95 in the HDF-S. The redshift and 
total magnitude information was obtained from two photo- 
metric redshift catalogs: Fernandez-Soto, Lanzetta & Yahil, 
(1999; HDF-N) and Labbe et al. (2003; HDF-S). The HDF- 
N catalog provides aperture magnitudes which are then cor- 
rected for the effect of neighboring. These magnitudes cor- 
respond to the SExtractor magnitudes named as "best" . For 
the HDF-S catalog, Labbe et al. provide with a "color" aper- 
ture (see their sec. 5.2). This "color" aperture is an isopho- 
tal aperture determined by the K s -band detection isophote 
at the 5 a detection threshold when the galaxies satisfy a 
"size" criteria, and is a fixed aperture when the galaxies are 
outside this "size" range. This "color" magnitude is finally 
transformed to a total magnitude using the total magnitude 
in the K 3 -band. The total K s -band magnitude is based on 
a Kron-like aperture (SEx AUTO) derived from the K s - 
detection image. When available, spectroscopic redshifts are 
used. This corresponds to 70% of the galaxies in the HDF-N 
and 60% in the HDF-S in the present work. 

The number identifications of the objects correspond to 
the same numbers used in the photometric redshift catalogs. 



3 THE FITTING ALGORITHM 

Our morphological analysis is based on the decomposition 
of the surface brightness (SB) profile into a bulge and a disc 
component. The bulge component is modeled by a Sersic 
(1968) profile of the form: 

I(r) = I b (0)e- b " ir/r ^ )1/n (1) 

where lb(0) is the bulge central intensity, r e ,b is the bulge 
semi-major effective radius and n the Sersic shape parame- 
ter. The intensity at the bulge effective radius r ej (, is given 
by Ie,b=Ib(0)e~ The quantity b n is a function of the 
shape paramenter n, and is chosen so that r Si b encloses 
half of the total luminosity. The exact value is derived from 
T(2n) = 27(2n,fe n ), where F(a) and j(a,x) are the gamma 
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and the incomplete gamma function (Abramowitz & Stegun 
1964, p. 260). The Sersic model contains the classical de Vau- 
couleurs profile when n is equal to 4. The disc component is 
an exponential profile (Freeman 1970) given by: 



J(r) = l d (0)e~ r/h 



(2) 



where ld(0) is the central intensity and h is the disc semi- 
major scale-length. The set of free parameters is completed 
with the ellipticities of the bulge e b and the disc e d - In addi- 
tion, it is possible to determine from the above parameters 
the total magnitude: 



hu{model) = -2.51og((F(, + F d )/t) + C 



(3) 



where t is the exposure time, C is the zeropoint, Ft is the 
bulge luminosity F b = I b (Q)r^ b (l - e b )2mT(2n) and 
Fd is the disc luminosity Fd = ld(0)h 2 (l — £d)27r. We have 
used C=23.02 (F606W) and C=22.09 (F814W) in the AB 
system. 

We can also evaluate the bulge-to-total luminosity ratio 
B/T: 



B 
T 



(4) 



Fb + Fd 

and the global semi-major effective radius r e by solving the 
equation: 



B/T 
1 - B/T 



2 1 {2n,b n {r e /r^ b ) L ' n ) 



T(2n) 



2 7 (2,r e /h)-l (5) 



Our parameter-recovering method is described in detail 
in Trujillo et al. (2001a) and Aguerri & Trujillo (2002). The 
main points of the routine are as follows. The final SB dis- 
tributions resulting from the convolution between the Point 
Spread Function (PSF) and our 2D Sersic model SB dis- 
tribution are dependent on the intrinsic ellipticity of the 
original source. However, to evaluate the intrinsic ellipticity 
of a model (i.e. the value of the ellipticity of the isophotes 
before being affected by the seeing) it is often insufficient to 
simply measure the ellipticity at one given radial distance. 
The ellipticity of the isophotes is reduced by seeing and this 
reduction depends on the radial distance of the isophote to 
the centre of the model, the size of the seeing, and the val- 
ues of the model parameters (Trujillo et al. 2001b, c). For 
that reason, the determination of the intrinsic ellipticity of 
the source and the fitting process to determine the struc- 
tural parameters must be done in tandem (i.e. iteratively 
and self-consistently) . To do this we simultaneously fit both 
the observed ID SB and ellipticity semi-major axis radial 
profiles using convolved profiles for each. We illustrate our 
technique in Fig. Q 

The measured SB and ellipticity semi-major radial pro- 
files are obtained from fitting ellipses to the isophotes of the 
selected objects. The elliptical isophotes are measured using 
the task ELLIPSE within IRAF 2 and fitting down to 1.5 a 
(where a is the standard deviation of the sky background 
of the images). A Levenberg-Marquardt non-linear fitting 
algorithm (Press et al. 1992) was used to determine the set 



2 IRAF is distributed by the National Optical Astronomical Ob- 
servatories, which are operated by AURA, Inc. under contract to 
the NSF. 



of free parameters which minimizes % 2 - A variety of start- 
ing parameters are used to ensure that our fits do not get 
trapped in local % 2 minima. 

The PSF was estimated fitting a Moffat function to star 
profiles. The high sampling of the drizzled images that we 
used and the presence of enough non-saturated stars along 
the field of view facilitates the estimation of the PSF prop- 
erties by a direct fit to this simple analytical function. We 
find that a Moffat function with the parameters presented 
below reproduces the observed PSF very well. We have es- 
timated a ~5 per cent uncertainty in the estimation of the 
FWHM due to changes from one WFPC2 position to an- 
other. In addition, the use of an analytical PSF facilitates a 
robust analytical estimation of the convolution between the 
PSF and the Sersic profile. 

We find the following /3 and FWHM (Full at Width 
Half Maximum) as our best-fitting stellar parameters: /9=3, 
FWHM=0'.'142 (F606W); 0=2.5, FWHM=0'.'147 (F814W). 
The best-fitting galaxy parameters in the F814W band are 
shown in Table Q There, we list the following: 

(1) The galaxy identification number as appearing in the 
photometric redshift catalogs Fernandez-Soto, Lanzetta & 
Yahil, (1999; HDF-N) and Labbe et al. (2003; HDF-S). 

(2) and (3) R.A. and Declination (J2000 coordinates). 

(4) Best AB(8140) model-independent magnitudes ac- 
cording to the photometric catalogs. 

(5) The surface brightness of the bulge at the effective 
radius (mag/arcsec 2 in the Isi4 band). 

(6) The semi-major effective radius of the bulge (arcsec) . 

(7) The Sersic index n of the bulge. 

(8) The ellipticity of the bulge. 

(9) The central surface brightness of the disc 
(mag/arcsec 2 in the Isi4 band). 

(10) The semi-major scale length of the disc (arcsec). 

(11) The ellipticity of the disc. 

(12) The bulge-to-total luminosity ratio (B/T < in- 
dicates galaxies classified as irregular/peculiar/merger sys- 
tems). 

(13) Spectroscopic redshift unless the index "a" appears 
which means photometric redshift from Fernandez-Soto, 
Lanzetta & Yahil, (1999; HDF-N) and Labbe et al. (2003; 
HDF-S). 

(14) Total AB(8140) magnitude as derived from the 
model fit. 

(15) The global semi-major effective radius of the galaxy 
(arcsec). 

For galaxies with z<0.4 the best-fitting galaxy param- 
eters in the F606W band are shown in Table [5] 

Our code has already been used and tested in a vari- 
ety of HST and ground-based studies of nearby and distant 
galaxies: the quantitative morphological analysis of galaxies 
in Abell 2443 (Trujillo et al. 2001a) and Coma (Gutierrez 
et al. 2004; Aguerri et al. 2004) clusters, the number den- 
sity evolution of galaxies in the SSA13 and SSA22 Hawaiian 
Deep Fields (Aguerri & Trujillo 2002), the ISAAC near- 
infrared images of the HDF-S (Trujillo et al. 2004) and the 
Flanking Fields HST images (Graham et al. 2004). In Tru- 
jillo et al. (2004) we have also made a comparison between 
our code and the GALFIT code (Peng et al. 2002) finding 
an excellent agreement. 
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Table 1. Best Morphological Parameter Values For Galaxies with Igi4(AB)<24.5 mag and z<l in the HDF-N and HDF-S. 



(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) (14) (15) 

ID R.A.(J2000) D.(J2000) I 8 i4(cat) /i e . b r e ,& n e b /x d (0) h e d B/T z I 8 i 4 (mod) r e 

(h:m:s) (o:m:s) mag mag/" 2 ("') mag/" 2 (") mag (") 



HDF-N 



1 

6 

7 

10 

11 

12 

18 

33 

54 

55 

81 

85 

91 

122 

124 

125 

138 

144 

183 

191 

197 

205 

218 

220 

245 

247 

257 

270 

273 

285 

310 

319 

335 

339 

344 

345 

350 

351 

363 

381 

391 

394 

409 

424 

445 

450 

454 

458 

465 

467 

476 

477 

495 

500 

512 

521 

524 

537 

549 

563 



12:36:57.47 
12:36:57.45 
12:36:57.60 
12:36:51.95 
12:36:56.42 
12:36:44.74 
12:36:49.24 
12:36:59.37 
12:36:58.65 
12:36:41.97 
12:36:45.30 
12:36:41.61 
12:36:42.28 
12:36:43.80 
12:36:46.50 
12:37:00.56 
12:36:52.01 
12:36:57.21 
12:36:41.33 
12:36:45.33 
12:36:45.05 
12:36:47.50 
12:36:46.26 
12:36:52.68 
12:36:46.54 
12:36:45.95 
12:36:51.72 
12:36:44.82 
12:36:50.16 
12:36:48.65 
12:36:48.12 
12:36:49.52 
12:36:53.45 
12:36:49.05 
12:36:50.85 
12:36:56.65 
12:36:58.77 
12:36:58.32 
12:36:49.99 
12:36:55.57 
12:36:41.94 
12:36:58.35 
12:36:55.62 
12:36:40.95 
12:36:58.07 
12:36:55.26 
12:36:56.93 
12:36:57.31 
12:36:43.62 
12:36:42.91 
12:36:47.28 
12:36:50.21 
12:36:41.48 
12:36:53.90 
12:36:54.78 
12:36:39.43 
12:36:50.26 
12:36:47.04 
12:36:55.16 
12:36:50.84 



62:12:10.5 
62:12:11.8 
62:12:12.3 
62:11:55.3 
62:12:09.2 
62:11:33.4 
62:11:48.4 
62:12:21.6 
62:12:21.7 
62:11:30.3 
62:11:42.8 
62:11:31.8 
62:11:34.7 
62:11:42.8 
62:11:51.3 
62:12:34.6 
62:12:09.6 
62:12:25.8 
62:11:41.0 
62:11:54.4 
62:11:54.1 
62:12:02.5 
62:11:59.7 
62:12:19.6 
62:12:03.1 
62:12:01.3 
62:12:20.1 
62:12:00.2 
62:12:16.9 
62:12:14.1 
62:12:14.8 
62:12:20.0 
62:12:34.2 
62:12:21.2 
62:12:27.2 
62:12:45.3 
62:12:52.3 
62:12:51.1 
62:12:26.2 
62:12:45.4 
62:12:05.4 
62:12:56.2 
62:12:49.2 
62:12:05.3 
62:13:00.3 
62:12:52.5 
62:12:58.1 
62:12:59.6 
62:12:18.2 
62:12:16.3 
62:12:30.6 
62:12:39.7 
62:12:14.9 
62:12:53.9 
62:12:58.2 
62:12:11.7 
62:12:45.7 
62:12:36.8 
62:13:03.6 
62:12:51.4 



21.2 
23.2 
23.1 
23.7 
23.2 
24.4 
22.9 
24.0 
23.3 
23.3 
23.9 
19.7 
23.8 
20.8 
21.9 
21.3 
22.9 
22.3 
22.4 
23.0 
24.2 
24.4 
24.4 
23.2 
24.3 
23.7 
21.5 
22.9 
22.4 
24.4 
24.4 
24.4 
22.8 
22.4 
24.2 
20.0 
21.2 
24.0 
24.4 
21.9 
21.0 
24.2 
22.9 
23.2 
22.3 
24.2 
23.5 
21.3 
22.7 
20.7 
22.7 
20.6 
23.4 
20.8 
24.3 
23.8 
21.4 
20.9 
23.5 
23.1 



20.3 
21.3 
24.2 
23.9 
0.00 
23.5 
0.00 
25.1 
0.00 
0.00 
23.5 
21.8 
22.8 
23.6 
20.9 
22.0 
21.8 
22.5 
23.7 
23.4 
26.4 
24.1 
0.00 
21.6 
26.1 
23.0 
0.00 
23.0 
0.00 
24.4 
22.2 
23.7 
24.2 
26.0 
24.2 
21.2 
21.9 
24.9 
0.00 
20.8 
0.00 
0.00 
23.9 
22.0 
21.8 
22.9 
0.00 
23.9 
0.00 
21.8 
22.8 
25.0 
20.1 
20.5 
24.1 
23.0 
22.2 
20.8 
21.7 
26.5 



0.10 
0.06 
0.08 
0.12 
0.00 
0.16 
0.00 
0.01 
0.00 
0.00 
0.53 
0.34 
0.15 
0.91 
0.16 
0.35 
0.14 
0.09 
0.09 
0.30 
0.12 
0.40 
0.00 
0.20 
0.55 
0.41 
0.00 
0.09 
0.00 
0.51 
0.13 
0.35 
0.25 
0.32 
0.39 
0.27 
0.08 
0.15 
0.00 
0.06 
0.00 
0.00 
0.12 
0.14 
0.25 
0.01 
0.00 
0.48 
0.00 
0.52 
0.38 
0.20 
0.06 
0.09 
0.34 
0.20 
0.23 
0.19 
0.14 
0.16 



2.91 
1.09 
1.38 
0.55 
0.00 
2.95 
0.00 
10.6 
0.00 
0.00 
0.30 
0.46 
0.27 
6.32 
2.30 
2.73 
0.53 
0.84 
10.7 
4.64 
8.75 
0.43 
0.00 
0.90 
0.68 
0.47 
0.00 
0.85 
0.00 
0.72 
0.77 
0.93 
1.28 
1.62 
0.49 
1.06 
1.46 
1.55 
0.00 
0.50 
0.00 
0.00 
0.61 
2.48 
0.89 
10.2 
0.00 
2.37 
0.00 
0.45 
1.08 
10.0 
1.44 
1.15 
0.51 
3.41 
6.31 
0.98 
2.05 
5.17 



0.20 
0.21 
0.71 
0.35 
0.00 
0.11 
0.00 
0.59 
0.00 
0.00 
0.64 
0.40 
0.64 
0.25 
0.13 
0.07 
0.28 
0.66 
0.10 
0.21 
0.23 
0.49 
0.00 
0.33 
0.38 
0.54 
0.00 
0.42 
0.00 
0.58 
0.35 
0.58 
0.78 
0.43 
0.33 
0.10 
0.64 
0.12 
0.00 
0.11 
0.00 
0.00 
0.46 
0.14 
0.10 
0.99 
0.00 
0.88 
0.00 
0.36 
0.68 
0.14 
0.69 
0.58 
0.32 
0.25 
0.22 
0.36 
0.29 
0.28 



20.99 
22.11 
22.91 
23.45 
21.70 
0.000 
20.92 
22.85 
20.69 
23.15 
0.000 
20.72 
23.26 
0.000 
0.000 
0.000 
23.42 
22.28 
21.24 
0.000 
23.21 
0.000 
23.28 
0.000 
20.74 
0.000 
19.77 
22.67 
20.86 
0.000 
0.000 
0.000 
20.94 
21.77 
0.000 
22.15 
21.36 
23.51 
21.74 
21.31 
19.90 
22.25 
23.56 
0.000 
23.78 
22.93 
21.25 
20.66 
20.70 
21.48 
23.51 
20.88 
20.66 
20.83 
0.000 
0.000 
21.46 
21.36 
0.000 
23.54 



0.36 
0.24 
0.36 
0.52 
0.33 
0.00 
0.21 
0.31 
0.14 
0.55 
0.00 
0.81 
0.34 
0.00 
0.00 
0.00 
0.61 
0.49 
0.26 
0.00 
0.30 
0.00 
0.39 
0.00 
0.09 
0.00 
0.30 
0.42 
0.41 
0.00 
0.00 
0.00 
0.31 
0.31 
0.00 
1.89 
0.48 
0.49 
0.23 
0.31 
0.26 
0.19 
0.62 
0.00 
0.95 
0.41 
0.21 
0.27 
0.19 
0.40 
0.68 
0.60 
0.13 
0.43 
0.00 
0.00 
0.38 
0.40 
0.00 
0.80 



0.34 
0.14 
0.26 
0.38 
0.61 
0.00 
0.37 
0.38 
0.30 
0.67 
0.00 
0.70 
0.26 
0.00 
0.00 
0.00 
0.46 
0.18 
0.39 
0.00 
0.07 
0.00 
0.65 
0.00 
0.50 
0.00 
0.63 
0.24 
0.80 
0.00 
0.00 
0.00 
0.62 
0.14 
0.00 
0.63 
0.32 
0.52 
0.78 
0.22 
0.06 
0.49 
0.31 
0.00 
0.60 
0.77 
0.54 
0.08 
0.31 
0.38 
0.25 
0.26 
0.43 
0.27 
0.00 
0.00 
0.38 
0.09 
0.00 
0.26 



0.36 
0.20 
0.01 
0.05 
0.00 
1.00 
0.00 
0.00 
0.00 
0.00 
-1.00 
0.14 
0.14 
1.00 
1.00 
1.00 
0.31 
0.02 
0.15 
1.00 
0.03 
-1.00 
0.00 
1.00 
0.32 
-1.00 
0.00 
0.05 
0.00 
-1.00 
1.00 
-1.00 
0.03 
-0.03 
-1.00 
0.19 
0.02 
0.10 
0.00 
0.09 
0.00 
0.00 
0.03 
1.00 
0.63 
0.00 
0.00 
0.05 
0.00 
0.62 
0.32 
-0.03 
0.28 
0.06 
-1.00 
1.00 
0.50 
0.32 
1.00 
0.01 



0.665 
0.561 

0.920° 

0.080" 
0.321 

0.880 a 
0.961 
0.472 
0.681 

0.200 a 
0.558 
0.089 

0.760 a 
0.764 
0.504 
0.562 
0.456 
0.561 
0.585 
0.372 

0.920 a 

0.920 a 

0.440 a 
0.401 
0.454 
0.679 
0.401 
0.457 
0.905 

0.800° 
0.962 
0.961 
0.559 
0.953 

0.640 a 
0.517 
0.320 

0.240 a 

0.010" 
0.790 
0.432 

0.520° 
0.950 
0.882 
0.319 

0.720 Q 
0.520 
0.473 
0.750 
0.454 
0.421 
0.474 

0.880 a 
0.642 
0.851 

1.000° 
0.678 
0.321 
0.952 
0.485 



21.1 
23.1 
23.4 
23.3 
23.1 
24.4 
22.8 
23.9 
23.3 
23.8 
23.7 
20.3 
23.7 
20.4 
21.9 
21.2 
22.7 
22.0 
22.5 
22.8 
23.8 
24.4 
24.4 
22.9 
24.2 
23.4 
21.4 
22.7 
22.5 
24.3 
24.5 
24.2 
22.4 
22.4 
24.2 
19.6 
21.3 
23.7 
24.5 
22.0 
20.8 
24.5 
22.9 
23.3 
21.8 
24.4 
23.4 
21.5 
22.6 
20.9 
22.2 
20.2 
23.3 
20.9 
24.4 
23.5 
21.3 
21.0 
23.3 
22.3 



0.38 
0.31 
0.59 
0.82 
0.50 
0.16 
0.37 
0.51 
0.22 
0.69 
0.53 
1.15 
0.48 
0.91 
0.16 
0.35 
0.62 
0.80 
0.40 
0.30 
0.49 
0.40 
0.59 
0.20 
0.21 
0.41 
0.46 
0.67 
0.58 
0.51 
0.13 
0.35 
0.51 
0.51 
0.39 
2.49 
0.79 
0.74 
0.38 
0.47 
0.46 
0.31 
1.01 
0.13 
0.39 
0.68 
0.35 
0.45 
0.31 
0.56 
0.80 
0.99 
0.15 
0.68 
0.34 
0.19 
0.46 
0.44 
0.13 
1.33 



HDF-N & HDF-S Quantitative Morphological Analysis 5 



(i) 

ID 



(2) 

R.A.(J2000) 



(3) 
D.(J2000) 



(4) 
I8l4(cat) 





(h:m:s) 


(' 


D:m:s) 


mag 


565 


12:36:48.99 


62: 


12:45.8 


23.6 


581 


12:36:49.47 


62: 


12:48.7 


23.9 


599 


12:36:57.71 


62: 


13:15.1 


23.0 


611 


12:36:50.82 


62: 


12:55.8 


22.5 


617 


12:36:38.97 


62: 


12:19.7 


22.2 


619 


12:36:55.46 


62: 


13:11.1 


22.1 


631 


12:36:55.15 


62: 


13:11.3 


23.6 


637 


12:36:51.43 


62: 


13:00.6 


22.7 


650 


12:36:44.18 


62: 


.12:40.3 


23.3 


653 


12:36:46.12 


62: 


12:46.4 


22.3 


655 


12:36:49.63 


62: 


12:57.5 


21.9 


662 


12:36:55.00 


62: 


13:14.7 


23.8 


669 


12:36:39.77 


62: 


12:28.5 


23.7 


671 


12:36:47.54 


62: 


12:52.6 


23.9 


683 


12:36:41.61 


62: 


12:35.6 


24.4 


689 


12:36:39.86 


62: 


12:31.5 


24.1 


694 


12:36:43.15 


62: 


.12:42.2 


21.7 


709 


12:36:44.18 


62: 


12:47.8 


21.6 


716 


12:36:38.40 


62: 


12:31.3 


22.6 


719 


12:36:49.68 


62: 


13:07.3 


24.2 


720 


12:36:43.96 


62: 


.12:50.0 


20.9 


727 


12:36:38.60 


62: 


12:33.8 


23.9 


740 


12:36:49.03 


62: 


13:09.7 


24.0 


743 


12:36:49.37 


62: 


13:11.2 


22.1 


746 


12:36:49.72 


62: 


13:13.0 


21.4 


749 


12:36:50.47 


62: 


13:16.1 


22.6 


751 


12:36:48.08 


62: 


13:09.0 


20.4 


757 


12:36:51.05 


62: 


13:20.6 


19.6 


774 


12:36:44.58 


62: 


.13:04.6 


21.2 


775 


12:36:48.47 


62: 


13:16.6 


23.2 


780 


12:36:48.78 


62: 


13:18.4 


22.6 


781 


12:36:46.75 


62: 


13:12.3 


23.6 


782 


12:36:45.61 


62: 


13:08.9 


23.6 


793 


12:36:49.07 


62: 


13:21.8 


24.4 


798 


12:36:46.15 


62: 


13:13.8 


24.2 


810 


12:36:42.51 


62: 


13:05.1 


24.3 


817 


12:36:42.72 


62: 


13:07.1 


22.1 


866 


12:36:47.45 


62: 


13:30.0 


23.5 


869 


12:36:45.86 


62: 


13:25.8 


20.9 


876 


12:36:45.41 


62: 


13:25.9 


22.6 


882 


12:36:52.23 


62: 


.13:48.0 


24.3 


884 


12:36:54.10 


62: 


13:54.3 


22.4 


888 


12:36:55.59 


62: 


13:59.8 


23.7 


890 


12:36:42.37 


62: 


13:19.3 


23.8 


899 


12:36:55.50 


62: 


.14:02.6 


22.9 


902 


12:36:42.02 


62: 


13:21.4 


23.5 


912 


12:36:51.80 


62: 


13:53.8 


21.1 


914 


12:36:49.44 


62: 


13:46.9 


18.1 


915 


12:36:48.03 


62: 


.13:43.0 


24.3 


938 


12:36:51.98 


62: 


.14:00.8 


23.0 


941 


12:36:52.87 


62: 


14:04.8 


23.2 


979 


12:36:53.65 


62: 


14:17.6 


23.4 


989 


12:36:49.50 


62: 


14:06.7 


21.8 


990 


12:36:46.03 


62: 


13:56.3 


24.2 


1018 


12:36:46.35 


62: 


14:04.6 


21.2 


1022 


12:36:51.39 


62: 


14:20.9 


23.3 


1027 


12:36:48.35 


62: 


14:12.3 


23.8 


1028 


12:36:50.36 


62: 


14:18.6 


23.0 


1029 


12:36:46.54 


62: 


14:07.5 


23.9 


1038 


12:36:48.10 


62: 


14:14.4 


24.4 


1048 


12:36:47.18 


62: 


14:14.2 


23.5 


1063 


12:36:53.11 


62: 


14:38.0 


24.2 


1067 


12:36:48.32 


62: 


14:26.2 


19.1 



(5) 
mag/" 2 



(6) 
r e ,b 
(") 



(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


n 


C6 


Md(0) 
mag/" 2 


h 

(") 




B/T 


z 


Isi4(mod) 
mag 


r e 
(") 


0.66 


0.48 


0.000 


0.00 


0.00 


-1.00 


0.512 


23.5 


0.41 


0.61 


0.03 


22.91 


0.21 


0.16 


0.45 


0.760" 


23.8 


0.19 


1.16 


0.05 


0.000 


0.00 


0.00 


-1.00 


0.952 


22.8 


0.51 


2.03 


0.37 


20.92 


0.30 


0.49 


0.13 


0.319 


22.1 


0.44 


0.98 


0.05 


0.000 


0.00 


0.00 


-1.00 


0.609 


22.2 


0.48 


3.98 


0.05 


0.000 


0.00 


0.00 


1.00 


0.370 


21.7 


0.72 


0.51 


0.17 


0.000 


0.00 


0.00 


-1.00 


0.321 


23.4 


0.23 


0.70 


0.23 


0.000 


0.00 


0.00 


-1.00 


0.089 


22.9 


0.89 


1.08 


0.11 


23.12 


0.47 


0.05 


0.03 


0.873 


22.7 


0.76 


2.20 


0.08 


22.73 


0.43 


0.25 


0.45 


0.900 


22.2 


0.33 


0.00 


0.00 


20.06 


0.26 


0.59 


0.00 


0.477 


21.9 


0.42 


0.37 


0.23 


22.02 


0.22 


0.63 


0.35 


0.511 


23.9 


0.29 


0.00 


0.00 


22.31 


0.26 


0.25 


0.00 


0.010 a 


23.6 


0.37 


0.00 


0.00 


22.07 


0.32 


0.66 


0.00 


0.681 


23.6 


0.55 


1.76 


0.24 


0.000 


0.00 


0.00 


1.00 


0.520™ 


24.3 


0.12 


0.66 


0.60 


23.07 


0.32 


0.52 


0.30 


0.640™ 


23.9 


0.34 


5.70 


0.25 


0.000 


0.00 


0.00 


1.00 


0.847 


21.2 


0.56 


1.79 


0.20 


21.25 


0.36 


0.20 


0.11 


0.558 


21.5 


0.53 


0.00 


0.00 


21.87 


0.58 


0.77 


0.00 


0.944 


22.6 


0.90 


0.85 


0.32 


22.51 


0.17 


0.18 


0.25 


0.920° 


24.2 


0.25 


1.79 


0.60 


20.58 


0.53 


0.60 


-0.06 


0.557 


20.8 


0.85 


2.85 


0.61 


21.95 


0.30 


0.66 


0.02 


0.904 


23.7 


0.49 


0.46 


0.44 


0.000 


0.00 


0.00 


1.00 


0.960™ 


23.8 


0.18 


2.50 


0.15 


0.000 


0.00 


0.00 


1.00 


0.478 


21.8 


0.19 


1.34 


0.70 


21.09 


0.57 


0.70 


0.13 


0.475 


21.4 


0.93 


17.0 


0.62 


22.39 


0.67 


0.71 


0.03 


0.851 


22.6 


1.09 


1.52 


0.04 


21.58 


0.57 


0.08 


0.41 


0.475 


20.2 


0.42 


0.58 


0.60 


20.99 


1.29 


0.60 


0.03 


0.199 


19.3 


2.10 


2.32 


0.72 


20.53 


0.53 


0.71 


0.05 


0.485 


21.1 


0.84 


0.99 


0.48 


23.15 


0.32 


0.09 


0.41 


0.474 


23.1 


0.32 


0.62 


0.14 


23.12 


0.64 


0.25 


0.02 


0.749 


22.3 


1.05 


0.77 


0.40 


22.32 


0.38 


0.63 


0.07 


0.600° 


23.4 


0.59 


2.20 


0.18 


0.000 


0.00 


0.00 


1.00 


0.480™ 


23.4 


0.17 


0.40 


0.05 


0.000 


0.00 


0.00 


-1.00 


0.880 a 


23.8 


0.30 


1.14 


0.31 


0.000 


0.00 


0.00 


1.00 


0.640™ 


24.1 


0.11 


3.00 


0.20 


0.000 


0.00 


0.00 


-1.00 


0.720 a 


24.5 


0.25 


10.0 


0.67 


19.82 


0.24 


0.74 


0.20 


0.485 


22.2 


0.44 


0.00 


0.00 


22.57 


0.26 


0.02 


0.00 


0.960 a 


23.5 


0.37 


1.55 


0.51 


21.96 


0.84 


0.28 


0.01 


0.321 


20.6 


1.39 


3.04 


0.67 


20.34 


0.18 


0.52 


0.43 


0.441 


22.2 


0.21 


0.68 


0.23 


0.000 


0.00 


0.00 


1.00 


0.520™ 


24.2 


0.19 


0.00 


0.00 


21.46 


0.32 


0.29 


0.00 


0.849 


22.3 


0.54 


6.45 


0.10 


22.41 


0.40 


0.72 


0.07 


0.559 


23.6 


0.62 


2.51 


0.24 


0.000 


0.00 


0.00 


1.00 


0.760 a 


23.8 


0.09 


0.00 


0.00 


21.48 


0.47 


0.80 


0.00 


0.564 


22.8 


0.80 


1.25 


0.64 


22.77 


0.20 


0.05 


0.46 


0.846 


23.6 


0.18 


0.58 


0.60 


22.25 


0.86 


0.60 


-0.49 


0.557 


20.8 


0.82 


2.40 


0.40 


0.000 


0.00 


0.00 


1.00 


0.089 


18.1 


0.55 


0.45 


0.60 


23.80 


0.55 


0.60 


0.04 


0.840° 


24.0 


0.88 


0.00 


0.00 


21.02 


0.32 


0.67 


0.00 


0.557 


22.7 


0.53 


1.52 


0.14 


21.96 


0.21 


0.05 


-0.08 


0.498 


23.3 


0.33 


0.00 


0.00 


21.36 


0.33 


0.73 


0.00 


0.517 


23.2 


0.48 


U.Dl 


AS 


91 9Q 

z± .zy 


U.OD 


n ^i 


u.zo 


n 7^9 

U. 1 oz 


91 f\ 
Zl .D 


n aa 


0.00 


0.00 


21.93 


0.16 


0.31 


0.00 


1.000° 


24.2 


0.25 


11.0 


0.05 


0.000 


0.00 


0.00 


1.00 


0.960 


20.6 


1.12 


0.00 


0.00 


22.05 


0.53 


0.60 


0.00 


0.439 


22.8 


0.88 


0.43 


0.80 


0.000 


0.00 


0.00 


-1.00 


0.920° 


23.9 


0.88 


0.79 


0.84 


22.45 


0.39 


0.37 


0.05 


0.816 


22.9 


0.62 


1.54 


0.33 


0.000 


0.00 


0.00 


1.00 


0.130 


23.7 


0.19 


0.57 


0.24 


22.03 


0.15 


0.37 


0.29 


0.920° 


24.3 


0.16 


0.53 


0.49 


21.06 


0.18 


0.43 


0.11 


0.609 


23.2 


0.27 


0.85 


0.60 


0.000 


0.00 


0.00 


-1.00 


0.400° 


24.0 


0.37 


1.02 


0.20 


19.40 


0.48 


0.20 


0.15 


0.139 


19.0 


0.68 



23.4 
22.6 
24.1 
22.3 
23.2 
24.3 
22.4 
25.0 
23.6 
21.0 
0.00 
23.7 
0.00 
0.00 
22.6 
21.9 
23.2 
21.7 
0.00 
24.2 
23.7 
24.4 
21.9 
21.3 
24.5 
22.6 
19.1 
21.7 
19.7 
22.2 
23.2 
23.2 
22.6 
23.5 
21.6 
24.5 
28.4 
0.00 
24.8 
20.1 
22.8 
0.00 
22.3 
21.6 
0.00 
21.0 
21.8 
19.4 
24.8 
0.00 
25.1 
0.00 
21.1 
0.00 
24.8 
0.00 
24.2 
23.1 
22.7 
21.8 
23.1 
23.5 
19.8 



0.41 
0.12 
0.51 
0.12 
0.48 
0.72 
0.23 
0.89 
0.08 
0.10 
0.00 
0.22 
0.00 
0.00 
0.13 
0.11 
0.57 
0.10 
0.00 
0.18 
0.37 
0.06 
0.18 
0.20 
0.77 
0.04 
0.10 
0.28 
0.05 
0.17 
0.07 
0.09 
0.18 
0.30 
0.11 
0.25 
1.88 
0.00 
0.32 
0.10 
0.10 
0.00 
0.03 
0.10 
0.00 
0.09 
0.56 
0.56 
0.15 
0.00 
0.19 
0.00 
0.18 
0.00 
1.13 
0.00 
0.88 
0.18 
0.20 
0.06 
0.15 
0.37 
0.18 



/. Trujillo and J.A.L. Aguerri 



(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) (14) (15) 

ID R.A.(J2000) D.(J2000) Isi4(cat) Me.b i e ,b n e b Md(0) h e d B/T z I 8 i 4 (mod) r e 

(h:m:s) (o:m:s) mag mag/" 2 ("') mag/" 2 (") mag (") 



HDF-S 



29 


22 


:32:46.27 


-60:34:16.3 


30 


22 


:33:00.48 


-60:34:17.6 


41 


22 


:32:56.06 


-60:34:14.1 


49 


22 


:32:45.60 


-60:34:15.1 


50 


22 


:32:47.55 


-60:34:08.5 


53 


22 


;32:50.94 


-60:34:15.2 


64 


22 


:32:55.21 


-60:34:10.0 


65 


22 


:33:01.54 


-60:34:10.7 


75 


22 


.32:55.23 


-60:34:07.4 
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Figure 1. Surface brightness and ellipticity semi-major radial profiles fitting to some galaxies in our sample. The three top galaxies are 
well fitted with a single Sersic law. The three bottom galaxies are fitted using a Sersic bulge plus an exponential disc. 



A preliminary way of examining the goodness of our 
structural parameter determination is to compare the best 
fit model magnitude with the magnitude measured in a 
model independent way (Fig. The model independent 
magnitude is the catalogs' magnitude. As expected from its 
extrapolation to the infinity, the model magnitude is almost 
always equal or brighter than the model independent deter- 
mination. 



4 OBJECT MORPHOLOGICAL 
CLASSIFICATION 

Following previous works (e.g. MS98) our galaxy classifica- 
tion scheme relies on B/T (Simien & de Vaucouleurs f986). 
Galaxies are considered "ellipticals" if their B/T>0.6. If this 
is the situation the galaxies are reanalysed with just a Sersic 



component and the parameters provided here correspond to 
this case. Galaxies with B/T between 0.5-0.6 are classified 
as SO and objects with B/T<0.5 are named "spirals". How- 
ever, there is a group of objects that can not be well fitted by 
our models. These objects have peculiar shapes (and conse- 
quently, their isophotes are hardly fitted by ellipses) and are 
classified in the Irr/Merger category. These objects are indi- 
cated in TablesQand|2]with B/T=-l. In addition, there are 
a group of visually classified spiral galaxies where a reliable 
fit is not achieved due to the presence of strong secondary 
structures as bars, spiral arms, etc. These galaxies are indi- 
cated with a value of B/T<0. 

The results of our galaxy classification are shown in 
Table 01 There are some discrepancies between the number 
of galaxies classified as E/S0 and Sa/Sc in the HDF-S and 
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759 


22:32:53.85 


-60:32:13.1 


24.2 


24.6 


0.30 


1.12 


0.87 


22.31 


0.27 


0.39 


0.06 


0.114 


23.6 


0.44 


763 


22:33:00.23 


-60:32:33.9 


20.5 


0.00 


0.00 


0.00 


0.00 


20.57 


0.45 


0.29 


0.00 


0.415 


20.7 


0.77 


766 


22:32:56.04 


-60:32:20.4 


22.5 


23.8 


0.21 


0.67 


0.23 


22.73 


0.72 


0.43 


0.05 


0.519 


21.9 


1.16 


768 


22:32:53.04 


-60:32:17.0 


23.7 


24.0 


0.21 


1.54 


0.29 


23.82 


1.04 


0.71 


0.16 


0.760° 


22.9 


1.45 


770 


22:32:48.87 


-60:32:16.0 


23.0 


25.1 


0.31 


2.08 


0.44 


21.41 


0.31 


0.61 


0.10 


0.580 


22.8 


0.51 


774 


22:32:52.01 


-60:32:14.9 


23.4 


0.00 


0.00 


0.00 


0.00 


22.00 


0.30 


0.39 


0.00 


0.800° 


23.2 


0.52 


822 


22:32:54.78 


-60:32:28.9 


24.1 


24.4 


0.11 


1.51 


0.10 


23.70 


0.41 


0.06 


0.08 


0.760° 


23.5 


0.64 


10010 


22:33:05.87 


-60:33:43.3 


22.6 


22.4 


0.37 


0.54 


0.27 


0.000 


0.00 


0.00 


-1.00 


0.640° 


22.4 


0.37 
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Figure 2. The total magnitude retrieved from the fitting model 
is compared to the total magnitude measured in a model indepen- 
dent way (catalogs' magnitude) for the galaxies analysed in this 
paper. Dashed lines corresponds to ±0.5 mag over the assumption 
of an exact agreement. 



Table 3. Coarse binning of galaxy classifications. 



Type 


B/T 


HDF-N 


HDF-S 


Total 


E/S0 


>0.5 


20% 


12% 


16% 


Sa/Sc 


<0.5 


63% 


74% 


68% 


Irr/Merger 




17% 


14% 


16% 



5 DETERMINATION OF ACCURACY 

To calibrate the accuracy of our parameter determination we 
have made two different tests. First, we have run our code 
on simulated artificial galaxies. This permits an evaluation 
of the apparent magnitude limit cut-off. Galaxies brighter 
than this magnitude limit are then claimed to be recovered 
with a certain degree of accuracy. We also have made a di- 
rect comparison between the structural parameters recov- 
ered using our algorithm and those recovered using GIM2D 
on the HDF-N (MS98). We describe these two issues in the 
following subsections. 



5.1 Simulations 



HDF-N. These differences illustrate the strong field-to-field 
variations present in small observing volumes. 



Simulations are a key element in the interpretation of the ob- 
servations. For fainter objects (i.e. at decreasing the signal- 
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Table 2. F606W Best Morphological Parameter Values For Galaxies with I 8 i4 (AB)<24.5 mag and z<0.4 in the HDF-N and 
HDF-S. 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


ID 


V606(cat) 




r e ,6 


n 


«6 


Md(0) 


h 


td 


B/T 


z 


V8i4(mod) 


r e 




mag 


mag/" 2 


(") 






mag/" 2 


(") 








mag 


(") 


HDF-N 


10 


24.4 


26.1 


0.11 


1.35 


0.26 


23.6 


0.44 


0.40 


0.02 


0.080™ 


23.9 


0.72 


11 


23.6 


0.00 


0.00 


0.00 


0.00 


22.1 


0.33 


0.60 


0.00 


0.321 


23.6 


0.50 


55 


23.9 


0.00 


0.00 


0.00 


0.00 


23.3 


0.53 


0.75 


0.00 


0.200° 


24.2 


0.75 


85 


20.1 


23.1 


0.28 


1.52 


0.10 


20.7 


0.70 


0.40 


0.06 


0.089 


20.0 


1.10 


191 


24.8 


25.5 


0.37 


4.40 


0.24 


0.00 


0.00 


0.00 


1.00 


0.372 


24.6 


0.37 


350 


21.8 


24.7 


0.15 


2.58 


0.55 


21.9 


0.53 


0.25 


0.02 


0.320 


21.6 


0.88 


351 


24.5 


27.3 


0.32 


2.10 


0.30 


23.7 


0.46 


0.54 


0.06 


0.240 a 


24.2 


0.74 


363 


24.9 


0.00 


0.00 


0.00 


0.00 


22.1 


0.18 


0.60 


0.00 


0.010 a 


24.8 


0.28 


445 


22.6 


21.7 


0.22 


0.87 


0.53 


23.4 


0.60 


0.56 


0.56 


0.319 


22.5 


0.37 


537 


21.5 


21.5 


0.24 


1.03 


0.42 


22.4 


0.52 


0.04 


0.36 


0.321 


21.4 


0.55 


611 


22.8 


0.00 


0.00 


0.00 


0.00 


20.5 


0.28 


0.70 


0.00 


0.319 


22.6 


0.48 


619 


24.0 


25.9 


0.72 


4.00 


0.23 


0.00 


0.00 


0.00 


1.00 


0.370 


23.9 


0.72 


631 


24.0 


22.6 


0.22 


0.37 


0.62 


23.4 


0.25 


0.25 


-0.52 


0.321 


23.9 


0.22 


637 


23.4 


26.3 


1.66 


1.00 


0.05 


0.00 


0.00 


0.00 


-1.00 


0.089 


22.6 


1.66 


669 


24.1 


0.00 


0.00 


0.00 


0.00 


22.6 


0.27 


0.29 


0.00 


0.010" 


23.9 


0.39 


757 


19.9 


21.7 


0.18 


0.42 


0.62 


21.3 


1.09 


0.51 


0.02 


0.199 


19.9 


1.79 


869 


21.5 


24.1 


0.23 


1.62 


0.67 


22.4 


0.64 


0.37 


0.03 


0.321 


21.8 


1.04 


914 


18.7 


20.3 


0.63 


2.30 


0.23 


0.00 


0.00 


0.00 


1.00 


0.089 


18.5 


0.63 


1029 


24.3 


23.2 


0.21 


1.86 


0.40 


0.00 


0.00 


0.00 


1.00 


0.130 


24.1 


0.21 


1067 


19.5 


20.0 


0.11 


0.51 


0.18 


19.2 


0.42 


0.34 


0.06 


0.139 


19.5 


0.66 


HDF-S 


184 


20.8 


23.6 


0.19 


6.03 


0.32 


20.7 


0.41 


0.27 


-0.08 


0.340 


20.9 


0.66 


186 


20.8 


23.7 


0.77 


6.84 


0.04 


0.00 


0.00 


0.00 


1.00 


0.328 


20.6 


0.77 


244 


21.7 


0.00 


0.00 


0.00 


0.00 


21.0 


0.67 


0.80 


0.00 


0.173 


21.7 


1.00 


283 


24.4 


27.2 


0.24 


3.11 


0.12 


23.6 


0.59 


0.55 


0.04 


0.180" 


23.6 


0.96 


314 


23.5 


0.00 


0.00 


0.00 


0.00 


21.6 


0.18 


0.14 


0.00 


0.220" 


23.5 


0.28 


438 


25.1 


0.00 


0.00 


0.00 


0.00 


22.6 


0.23 


0.55 


0.00 


0.359 


24.7 


0.37 


446 


21.9 


24.0 


0.12 


12.1 


0.09 


21.2 


0.23 


0.44 


0.19 


0.364 


22.8 


0.37 


487 


24.2 


24.2 


0.43 


0.69 


0.14 


0.00 


0.00 


0.00 


-1.00 


0.240" 


23.6 


0.43 


508 


24.0 


24.3 


0.37 


0.64 


0.51 


0.00 


0.00 


0.00 


-1.00 


0.280" 


24.8 


0.37 


528 


21.7 


0.00 


0.00 


0.00 


0.00 


21.1 


0.45 


0.56 


0.00 


0.365 


21.8 


0.68 


676 


24.5 


23.4 


0.25 


1.42 


0.39 


0.00 


0.00 


0.00 


-1.00 


0.318 


24.1 


0.25 


759 


24.4 


23.8 


0.18 


0.36 


0.26 


23.0 


0.31 


0.48 


0.17 


0.114 


24.0 


0.42 



to-noise S/N) the structural parameters of the galaxies are 
progressively more difficult to determine. To quantify this ef- 
fect it is necessary to run realistic simulations matching the 
observational conditions. We simulated two kind of objects. 
First, we created bulge-only (i.e. pure elliptical) artificial 
galaxies, and second, we constructed bulge+disc galaxies. 
In both cases, we made 150 artificial objects with structural 
parameters randomly distributed in the following ranges: 

• bulge-only structures: 19<I(AB)<25, 0.05"s:r e <0.6", 
0.5s: n sC 4, and < e sC 0.6 

• bulge+disc structures: 19<I(AB)^25, 
0.05"s:re sCO.6", 0.5s: n < 4, and < e b < 0.4, 
0.2" sC h sC 1.5", Os: B/T sC 1, and sC e d sC 0.6. 

Artificial galaxies were created by using the IRAF task 
MKOBJECT. We support as an input to this task the sur- 
face brightness distribution coming from our detailed convo- 
lution between the PSF and the original model. To simulate 
the real conditions of our observations, we added a back- 
ground sky image (free of sources) taken from a piece of the 



real image. The PSF FWHM in the simulation was set at 
0.147" and assumed known exactly. The pixel scale of the 
simulation was 0.04", as are the drizzled HDF pixel scales. 
The same procedure was used to process both the simulated 
and actual data. 

From these simulations we find that the bulge and disc 
parameters (as well as the galaxy global properties as to- 
tal magnitude, r e , etc.) can be determined with an accu- 
racy of ~10% for galaxies with Igi4(AB)<24.5 mag. Con- 
sequently, we consider Isi4(AB)=24.5 mag as our limiting 
magnitude for a reliable determination of the parameters. 
Further details of how the simulations are constructed and 
the errors determined can be found in Trujillo et al. (2001a) 
and Aguerri & Trujillo (2002). 

5.2 Comparison with Previous Works 

MS98 using a different code, GIM2D (Simard 1998), con- 
ducted a quantitative morphological analysis on the HDF-N 
close to that presented here. The main differences between 
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both analyses are: the limiting apparent magnitude and the 
model for the bulge component. GIM2D was run down to 
i8i4(AB)<26 mag, i.e. 1.5 mag deeper than in the present 
analysis. However, although individual errors for the galaxy 
parameters are available, there are not simulations showing 
the accuracy of the structural parameters as a function of 
the apparent magnitude. Another difference is that GIM2D 
was run using a de Vaucouleurs model (n=4 fixed) for the 
bulges. Fixing n to 4 has the advantage of decreasing the 
number of free parameters, however, in the nearby universe 
a range of n (from 0.3 to 6) is found (Andredakis, Peletier & 
Balcells 1995, Graham 2001). Moreover, n relates with the 
total magnitude of the bulge component and it would be of 
interest to know if this relation evolves with time (Aguerri 
& Trujillo 2004). 

In Fig. [3] we show the comparison between the global 
size of the galaxies x e as obtained from our code (n free) and 
GIM2D (n=4) versus z and versus the apparent magnitude. 
The comparison is limited to i8i4(AB)<24.5 mag and z<l. 
There is a good agreement between these two different algo- 
rithms (~75% of the galaxies have a size difference less than 
20%). However, we find a small offset of 5% toward smaller 
r e in our code relative to GIM2D. This offset seems to be 
real in the sense that a T-Student test rejects the hypothesis 
that this is produced by chance by more than 99.9%. This 
offset is not dependent on z or magnitude and (due to it is 
not present in our simulations) must be intrinsically related 
to the way our code and GIM2D estimate the global galaxy 
size. As explained before our code uses Eq. (5) to estimate 
r e . This equation is derived by the definition of effective ra- 
dius (i.e. L(r e )=0.5L(oo)). On the other hand, GIM2D eval- 
uates r e by integrating the bulge and disc expressions out 
to infinity and then evaluating r e . It is possible that their 
assumption of n—4 plays here this minor role including more 
light than what it is actually observed with a n smaller. If 
this is the case we would expect that their size would be 
slightly larger than ours. In any case, it is important to note 
that this offset is much smaller than the standard deviation 
found at comparing both codes (~20%). 

We also present the comparison between our B/T ratio 
and GIM2D in Fig. 2] The global properties of the galaxies 
r e and B/T are shown to be very robust for this sample of 
galaxies. In the nearby Universe, Simien & de Vaucouleurs 
(1986), using a n=4 bulge, did an estimation of the B/T 
ratio for a sample of galaxies ranging a large variation of 
morphological types. Graham (2001) conducted a similar 
analysis using this time a n free bulge. He found that the 
use of a n=4 bulge (instead of a n free) produces an overes- 
timation of the B/T ratio. We do not see a strong difference 
between our B/T values and the GIM2D estimations. This 
is probably due to most of our galaxies having small B/T 
ratios and consequently the bulge component is not domi- 
nant. However, the use of a different model for the bulge will 
have consequences when dealing with bulge properties such 
as r e ,j,. 

We also compare our B/T determination with the mor- 
phological type obtained from van den Bergh et al. (1996, 
hereafter vdB96). This is shown in Fig. |K| The numerical 
classification used by vdB96 (Abraham et al. 1996) is as 
follows: E/star: -1, E: 0, E/SO/Sa: 1, SO/Sa: 2, Sa/Sab: 3, 
Sb/S/Ir: 4, Sc: 5, Ir: 6, Peculiar/Mergers: 7 or 8. As ex- 
plained before, galaxies with B/T>0.6 were reanalysed and 



Table 4. Coarse comparison of galaxy classifications 
in the HDF-N down to l8i4(AB)<24.5 mag and z<l 
for the same group of galaxies using three different 
methods. The number of analysed galaxies is 102. 



Type 


B/T 


This work 


MS98 


vdB96 


E/S0 


>0.5 


20% 


15% 


22% 


S/Irr 


<0.5 


80% 


85% 


78% 



1 .0 




d(r e )/r e 



Figure 3. The relative error between the size estima- 
tion using our code and GIM2D d(r e )/r e =2x(r e (GIM2D)- 
r e (ours)) / (r e (GIM2D)+r e (ours)) is shown versus z and versus the 
apparent Isi4 magnitude. The histogram shows that for ~75% of 
the galaxies the difference on size is less than 20%. 

fitted with a bulge-only Sersic component (i.e. B/T=l). For 
illustrative purpose, to avoid the appearance of just one 
point at B/T=l in Fig. |S] we have randomly moved all the 
galaxies with this B/T value a 5% around 0.95. The main 
difference between our analysis and MS98 is the group of 
galaxies classified by vdB96 as vdB=0. Contrary to MS98, 
most of the vdB96 visual ellipticals are also classified as 
B/T=l in our quantitative analysis. 

The fractions of the different morphological types ac- 
cording to vdB96 (down to Isi4(AB)<25 mag and not z 
restrictions) are: 30% E/S0, 31% S/Irr and 39% not clas- 
sified. However, these percentages depend on the limiting 
magnitude and the limiting redshift used. To make a reason- 
able comparison between the different morphological types 
which follow from the three different methods (ours, MS98 
and vdB96) we have used exactly the same galaxies. Our 
results are presented on Table [3] As we can see the agree- 
ment between the different methods is good. Consequently, 
previous claimed discrepancies between visually and quanti- 
tative methods (MS98) are shown here to be potentially due 
to different limiting magnitude criteria and redshift. How- 
ever, although the fractions may be similar for each study, 
different galaxies may populate each bin. We have checked 
if this is the case by comparing galaxy by galaxy and found 
that the agreement in each bin is greater than 65%. 

We have also undertaken a comparison with the anal- 
ysis done in the HDF-S by Moth & Elston (2002). In that 
work, 33 galaxies with 0.5^z^l.2 and Isi4<25 are analysed. 
The model used for extracting structural information is a 
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Figure 4. The difference between the B/T estimation using our 
code and GIM2D d(B/T)=B/T(GIM2D)-B/T(ours) is shown ver- 
sus z and versus the apparent L314 magnitude. The histogram 
shows that for ~73% of the galaxies B/T differs less than 0.2 
units. 




Figure 5. Comparison of the parameter B/T with the vdB mor- 
phological type derived from the visual classification of van den 
Bergh et al. (1996). 



single Sersic law and not a bulge+disc decomposition. For 
that reason the comparison is not straightfoward. For the 20 
galaxies we have in common, we find a very good agreement 
in the global size (difference less than 20%) for 18 objects 
(90%). A similar analysis using also a single Sersic law is 
done in Trujillo et al. (2004) for K s -band selected galaxies 
down to K IS =23.5. Trujillo et al. analysed in the Igi4-band 
galaxies up to z~0.8. The agreement in the global sizes for 
the common 61 galaxies is good, with ~68% of the objects 
having a difference less than 10% . 



6 COMPLETENESS AND 

REDSHIFT-DEPENDENT SELECTION 
EFFECTS 

The ability to detect a galaxy is critically dependent on its 
apparent surface brightness. Galaxies with larger B/T are 
easier to detect because they are more concentrated, and 
galaxies with larger radii are harder to detect than smaller 



ones at given total magnitude. Therefore, it is of interest to 
know whether our observed parameter distribution would 
be biased by the galaxy selection function of the SExtractor 
detection algorithm (used in both catalogs we are using). 
The evaluation of completeness of the HDF-N was carefully 
model in detail by MS98 (their Fig. 5). According to these 
authors galaxies brighter than Isi4(AB)<24.5 are fully re- 
covered. Due to the very similar characteristics of the HDF- 
S these results are also valid for that field. This means that 
in the sample of galaxies shown in this work the incomplete- 
ness is not playing any role and consequently we will not 
take it into account in the subsequent analysis. 

Redshift -dependent observational biases can also mimic 
real evolutionary changes in the galaxy population (Simard 
et al. 1999). For a given magnitude limit (in our case 
Igi4(AB)=24.5 mag) there is a corresponding threshold in 
the observable restframe luminosity which increases with 
redshift. We illustrate this in Fig. HJ This figure shows the 
rest-frame luminosity in the V-band versus the redshift. 

To estimate the rest-frame luminosity in the V-band 
we make use of the complete filter coverage of the HDFs. 
We have used the fluxes measured in the U300, B450, V606, 
I814, J, H and K bands provided by the Fernandez-Soto et 
al. (1999; HDF-N) and Labbe et al. (2003; HDF-S) cat- 
alogs. With these fluxes we construct observed SEDs for 
every galaxy. As we have redshift data for our sources we 
estimate the redshifted wavelength of the center of the local 
V-band filter (A=5510 A) (i.e. we evaluate A=5510x (1+z)). 
At that wavelength, we measure the restframe V-band flux 
by linearly interpolating the catalog's fluxes at the two clos- 
est A values. The advantange of using this method is that we 
do not rely on any SED model to estimate the luminosity. 
There are two main sources of uncertainty: 

(i) The uncertainty due to the photometric redshifts: 
most of our galaxies (70% in the HDF-N and 60% in the 
HDF-S) have z spectroscopically determined, however, the 
remaining fractions have an uncertainty of dz/(l+z)~0.1. 
This uncertainty in z can be translated into luminosity un- 
certainty of ~30%. 

(ii) The uncertainty due to the total flux calibration: the 
reliability of the total flux measured through a given fil- 
ter will depend on the S/N of the sources. Fortunately, the 
HDFs are so deep, and the objects selected in this work so 
bright, that this effect does not play a major role. In fact, we 
can see from the comparison between our model magnitudes 
(which are an extrapolation of the total magnitude up to an 
infinite radial extension) and the aperture magnitude (our 
Fig. |2l that the mean discrepancies between both is ~0.2 
mag, being the "true" flux value probably in between these 
two estimates. The above uncertainty in the magnitude is 
translated into a luminosity uncertainty of ~20%. 

According to the above two points we will assume on 
what follows an uncertainty in the restframe V-band lumi- 
nosities of ~35%. 

In Fig. [(J we show how our high-z sample represents 
only the most luminous fraction of the galaxy population. 
The lack of bright galaxies at low redshift results directly 
from their low volume densities coupled with the small vol- 
ume of the HDF over this redshift range. Galaxies with 
h v > 2 x 10 9 h" 2 L (i.e. M B <-17.5 assuming B-V=0.9) 
can be studied along the full redshift range. 
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Figure 6. The hy-z diagram for the objects selected in our sam- Figure 8. The intrinsic (i.e. the recovered non-PSF affected) 

pie with f8i4(AB)<24.5 mag. cllipticity of the galaxies versus: a) the redshift of the observed 

sources and b) the apparent Igi4 total magnitudes. No clear re- 
lation is observed. 
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Figure 7. The relation between the model-independent global 
magnitude versus the apparent global semi-major effective radius 
size (in the Igi4— band). 



7 HDF-N AND HDF-S LUMINOSITY-SIZE 
RELATIONS OUT TO Z~l 

It is of interest to know whether the population of galaxies 
in the HDF-N and HDF-S present similar global structural 
properties. We show in Fig. Q the relation between the 
model-independent global magnitude versus the apparent 
global semi-major effective radius size (in the Isi4-band). 
We have run the generalization of the Kolmogorov-Smirnov 
(K-S) test on two dimensional distributions (Fasano & 
Franceschini 1987) and found that the null hypothesis (i.e. 
that the two distribution have a similar origin) can not be 
rejected. 

In what follows, to facilitate the comparison with local 
calibrating data, we use, instead of the global semi-major 
effective radius r e , the circularized effective radius r e , c (i.e. 
we use r e , c = r e \l (1 — e) with e the intrinsic, non-PSF af- 
fected, projected ellipticity of the galaxy). For those galaxies 
where B/T<1, e — e<j. To be sure that the estimate of the in- 
trinsic ellipticity of our sources, and hence the conversion to 
circularized effective radius, was not systematically affected 
by the PSF effect we studied the presence of any obvious 



trend with z or the Isi4 apparent magnitude (Fig. |HJ . No 
clear relation is found. 

We have conducted a similar analysis to that presented 
in Fig.Qusing the rest-frame properties of the above galax- 
ies. This is shown in Fig. [5^,. In this case, the null hypoth- 
esis also can not be rejected. To avoid any bias due to K- 
corrections we have used for galaxies with z <0.4 the struc- 
tural parameters obtained in the F606W filter (Table 
For galaxies with 0.4< z <1 the filter used was F814W. 
Therefore the sizes are computed for roughly the same rest- 
frame wavelenghts (between B and V) over the entire red- 
shift range. This avoids any band-dependent size effects. 
This would be particularly important in late-type galaxies. 
In the local universe, Graham (2001) finds (assuming that 
the effective radius of the bulge is not modified) that the 
scale of the discs increases 15-25% from R to B band. 

In Fig. Eft we show the same sample but making 
a separation between early-type (B/T>0.5) and late-type 
(B/T<0.5) galaxies. At a fixed luminosity, late-type galax- 
ies are larger than the early-types. 

7.1 Testing the hypothesis of evolution 

Shen et al. (2003) have shown the median and dispersion of 
the distribution of the circularized half-light radius (Blan- 
ton et al. 2003) for different bands as a function of the lu- 
minosity in the Sloan Digital Sky Survey (SDSS, York et 
al. 2000) . Their analysis is based on an unprecedented large 
database of nearby galaxies with typical z~0.1. In the fol- 
lowing, we closely match the analysis presented in Trujillo 
et al. (2004). We use the SDSS g-band (the closest avail- 
able filter to our V-band) size distributions (S. Shen, 2003, 
private communication) as a local reference of the size dis- 
tribution of galaxies in the nearby universe. The median and 
dispersion of the Shen et al. distribution in the g-band sep- 
arating into early and late type are shown in Fig. 1101 We 
also overplot our HDFs data. 

The criteria followed to separate between early and late 
types in the presented SDSS relations is different from the 
one we have used. In the Shen et al. distributions a galaxy 
is separated into early or late type according to whether its 
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Figure 9. a) The distribution of circularized rest— frame sizes 
versus the rest-frame V— band luminosities is shown separating 
galaxies according to their ownership to the HDF-N or to the 
HDF-S. b) Same as a) but this time segregating the points be- 
tween late-type (B/T<0.5) and early— type (B/T>0.5) galaxies. 
For clarity error bars are not shown. The mean size relative error 
is 20% and the mean luminosity error 30%. 
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7.1.1 Late-Type galaxies 

At a given luminosity, Shen et al. propose for nearby galaxies 
the following size log-normal distribution with median f e 
and logarithmic dispersion a: 



f(r e \f e (L),a(L)) = 



1 



27rt(L) 



(6) 



To address the problem of size evolution, we have as- 
sumed that the functional form of the size log-normal distri- 
bution function (Eq. [SJ found for nearby galaxies also holds 
at all redshifts, but with evolving parameters: 



f e (L,z) =f e (L,0)(l + z)- a 
a{L,z)=a{L,0)(l + z) rJ . 



(7) 
(8) 



Here, f e (L,0) and a(L,0) are the median size and disper- 
sion provided at z=0 by the Shen et al. (2003) data, and a 
and (3 describe the redshift evolution. Implicit in Eq. [7|and 
Eq |5] we assume the same size evolution for all the galax- 
ies independently of their luminosity. However, in order to 
assure that we observe the same type of galaxies along the 
explored redshift range 0<z<l (see Fig. |SJ we select only 
galaxies with Ly>2xl0 9 h -2 Lq. To explore in much de- 
tail the evolution of the luminosity-size relation with z we 
have separated our sample in galaxies with z<0.5 (mean 
z~0.4) and galaxies with z>0.5 (mean z~0.7). The number 
of galaxies with Lv>2xl0 9 h -2 L© is 31 in the lower z bin 
and 85 in the higher. Due to the scarce number of objects 
in the lower z bin we will concentrate our analysis on the 
galaxies with z>0.5. In addition, there is a group of out- 
liers with z<0.5, luminosities between 0.1-lxlO 10 h -2 Lq 
and small sizes that affects our analysis at lower redshifts. 
These galaxies have all z photometrically determined. The 
option of taking only galaxies with z spectroscopically de- 
termined for this bin is not useful because of the few number 
of galaxies then available. 

We evaluate a and f3 using a Maximum Likelihood 
Method. The likelihood to maximize is given by: 



n/ 



1 9L <X> flr e 



0) 



Figure 10. a) The distribution of rest-frame sizes versus the 
rest— frame V— band luminosities is shown for late— type (B/T<0.5) 
galaxies, b) Same as a) with early— type (B/T>0.5) galaxies. The 
points are separated according to the z value. For clarity error 
bars are not shown. The mean size relative error is 20% and the 
mean luminosity error 30%. 



global Sersic index n is larger or smaller than 2.5. The index 
n derived for these authors corresponds to a global Sersic fit 
to the whole galaxy, whereas we have used a bulge+disc de- 
composition and our index n refers only to the bulge compo- 
nent. The separation between early and late type according 
to the global index n can be affected by the inclusion of 
some Sa galaxies as early-type galaxies. This effect is well 
illustrated by Ravindranath et al. (2004, their Fig. 1). Also, 
some contamination in the late-type galaxies at lower lumi- 
nosities is expected from dwarf ellipticals which have n~l. 



Each member (one per galaxy measured) of the product is 
a convolution between the expression / given in Eq. HJ and 
a gaussian representing the error on the luminosity ql and 
another gaussian representing the error on the size g Te ■ The 
width of each gaussian is given by the error at estimating 
the luminosity and the size respectively. 

The result of this analysis is shown in Fig. 1111 Ac- 
cording to this analysis the luminosity-size distribution of 
disc galaxies presents a slightly (although compatible with 
cero) broader dispersion at higher redshift. The dispersion 
increases by ~ 20(±25)%. To know whether this increase in 
the dispersion is real or not will require a large number of 
galaxies. 

We observe a moderate decrease in the sizes ~ 
30(±10)% at z>0.5. In the luminosity-size distribution a 
decrease in the size can also be interpreted as an increase 
in the luminosity of the objects. In fact, this seems to be a 
plausible hypothesis to explain the observed luminosity-size 
evolution found at high-z. According to the studies of the 
evolution of the stellar mass-to-light ratio (see e.g. Rudnick 
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Figure 11. Luminosity-size evolution for late— type galaxies with 
z>0.5. Likelihood contours representing the evolution in size and 
dispersion in the a— ft plane. Solid line represents the lcr con- 
fidence level contour, dashed line the 2a confidence level. The 
cross shows the position of no— evolution in this plane. Positive 
values of a represent decreasing values of the size with redshift. 
Positive values of /3 represent increasing the intrinsic dispersion 
a of the population with redshift. The luminosities are given in 
10 10 L n . 



Figure 12. The distribution of rest— frame sizes versus the rest- 
frame V— band luminosities is shown for early— type (B/T>0.5) 
galaxies. The points are separated according to the z value. Small 
open points correspond to 200 ellipticals from a nearby sample 
of Virgo, Fornax and Coma clusters. For clarity error bars are 
not shown. The mean size relative error is 20% and the mean 
luminosity error 30%. 



et al. 2003) the stellar populations were younger at high-z, 
and for that reason, the mass-to-light ratio is a decreasing 
function with z. Consequently, even though the size of the 
galaxies does not decrease with z, the luminosity evolution 
of the objects will simulate this effect. Consequently, our ob- 
served evolution in the luminosity-size distribution can be 
alternatively interpreted as a decrease in the central surface 
brightness of ~ 0.77(±0.30) mag at at z~0.7. 

The results presented here are in agreement with those 
found by Lilly et al. (1998), Simard et al. (1999) and Ravin- 
dranath et al. (2004). These authors do not find a great 
change from local galaxies out to z—0.7 and a hint of change 
at z>0.7 with shorter objects than those locally observed for 
a similar magnitude range. 

7.1.2 Early-Type galaxies 

In addition to the SDSS distribution and for covering the 
range of faint galaxies in the relation shown in Fig. 110b . 
we have used a local sample of 200 elliptical galaxies from 
the Virgo, Fornax and Coma cluster (Caon et al. 1990; 1994; 
Binggeli & Jerjen 1998; Gutierrez et al. 2004). This sample 
ranges from -22<M_b<-12. To transform to V-band magni- 
tudes we have used B — ^=0.9. To facilitate the comparison 
with this data, this time we use the semi-major effective 
radius. The result of this comparison is shown in Fig. 1121 

It is interesting to note that the observed distribution, 
with a plateau at the faintest galaxies, is in qualitative good 
agreement with the distribution of galaxies found in nearby 
samples. The apparent absence of evolution in the shape of 
the distribution of galaxies in the size-luminosity relation 
suggests that the early-type galaxies were yet assembled 
at higher redshift (Fasano et al. 1998, van Dokkum & Ellis 
2003, Franx et al. 2003, van Dokkum et al. 2003, Daddi et al. 
2003, Vazdekis, Trujillo & Yamada 2004). However, we ob- 
serve that at a given size, high-z early-type galaxies present 



higher luminosities than the local sample. To quantify this 
increase we have fitted the early-type luminosity-size rela- 
tion of the SDSS by using an expression of the form: 

r e = axL b v (10) 

We find the following values: a5_oss=2.25 and bsrjss=0.65. 
We have repeated the same exercise this time in our z>0.5 
early-type sample. Due to the scarcity of points we have 
fixed the b slope to the same value than the SDSS sam- 
ple. To make the fit we use the points defined by the gi- 
ant ellipticals (i.e. those points with Ly>10 10 h~ 2 L©). On 
doing that we find: aHOFs=l(±0.05). This implies a fac- 
tor of ~3.5(±0.25) greater in luminosity which corresponds 
to ~1.35(±0.10) mag. This is in good agreement with the 
brightening of 1-1.5 mag observed in the rest-frame B-band 
by Im et al (1996; 2002) and Schade et al. (1999) using the 
HST Medium Deep Survey and other HST fields. It is also 
obtained in large area fields (Bell et al. 2004). It is of in- 
terest to note that the comparison between SDSS and the 
HDFs are between field early-type galaxies. The bright el- 
liptical galaxies in the nearby clusters define a relation with 
SL c iusters~3A5. Consequently, we find that early-type galax- 
ies in the HDFs at z~0.7 are a factor of 6.7 brighter in 
luminosity (~2 mag) than the present-day elliptical in clus- 
ters. In agreement with previous studies that indicate that 
early-types in clusters are older than in the field. 

We have checked, using a 2D KS, test whether the 
no-evolution (versus the passive evolution) model can be 
ruled out. To do that we have used only the galaxies with 
Ly>0.1xl0 10 h~ 2 L and z>0.5. This test rejects the no- 
cvolution hypothesis with a confidence larger than than 
99.9%. This luminosity evolution is consistent with models 
in which the bulk of stars in the field E/S0 galaxies formed 
before z>1.5 and have been evolving rather quiescently. 
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8 DISCUSSION AND CONCLUSIONS 

It is of interest to compare the results presented in this paper 
with the theoretical expectations according to the different 
scenarios of galaxy formation. 

In the hierarchical scenario, following Mo et al. (1998), 
the size of the baryonic component of disc galaxies is ex- 
pected to scale with redshift as R oc H~ 1 (z) at fixed circular 
velocity, or R oc H~ 2 ^ 3 (z) at fixed mass. In a flat universe 
0=1, the Hubble constant evolves with redshift as: 

H(z) = H [n A + Q {l + z) 3 ] 1/2 (11) 

where Qa and Qo are the present-day density parameters. 
At z~0.7 (the mean z of our high-z population), the ex- 
pected evolution is 0.68<R(0.7)/R(0)<0.77. This decrease 
is compatible with our estimation: ~30(±10)%. However, 
before extracting any conclusion we must take into account 
that the Mo et al. predictions are for the mass-size relation 
whereas this work deals with the luminosity-size relation. In 
fact, according to previous (Trujillo et al 2004) and ongoing 
(Barden et al. 2004; private communication) work, although 
the luminosite-size relation has changed dramatically from 
z~3, the stellar mass-size relation has remained nearly in- 
variant from then. 

In the infall scenario, for the rest-frame B-band, 
Bouwens & Silk (2002) propose the following dependence 
in size: R(z)/R(0) = l-0.27z. At z—0.7 this corresponds to 
a decrease of 19%. This value is in agreement (although is 
slightly smaller than the observed) with our measurements. 
Contrary to the above predictions of the hierarchical mod- 
els, the Bouwens & Silk prediction is for the luminosity-size 
relation, and consequently, the comparison with our obser- 
vations is a more direct one. It is important to note, however, 
that our analysis is done in the V-band, and the expected 
evolution from the Bouwens & Silk (2002) model could be 
slightly different in this band. 

As pointed out by Bouwens & Silk (2002): "The hier- 
archical and infall models predict relatively similar amounts 
of evolution in global properties (size, surface brightness, 
and luminosity) for disc galaxies to 2~1". This implies we 
need to find an extra test for disentangle between both 
scenarios. The main difference between both approaches is 
the evolution of the comoving number density of the ob- 
jects. In the infall models is implicitly assumed no num- 
ber evolution whereas in the hierarchical model the number 
of disc galaxies, is an increasing function of z: N(V c ,z)dV c 
oc (1 + £)( 2 ' 9-3 ' 3 ) at a fixed circular velocity, or N(L,z)dL 
oc (1 + af )(i- 4 s-i-85) ) a t a fixed luminosity (see e.g. Mao, Mo 
& White 1998). The scarce number of galaxies in this work 
and the strong field-to-field variations prevent a definitive 
answer at this point and larger areas are required. 

Contrary to the disc galaxies, the number of field E/SOs 
decreases at high-z in the hierarchical picture. However, in a 
flat nonzero A universe the semianalytical models predict a 
weak number density evolution. This makes difficult to dis- 
entangle between the monolithic evolution scenario and the 
hierarchical view up to «~1. A more definitive answer to this 
problem will come from studying in detail the population of 
E/SOs galaxies in the l<z<2 interval. 



ACKNOWLEDGMENTS 

We are happy to thank Shiyin Shen for providing us with 
the Sloan Digital Sky Survey data used in this paper. We 
are grateful to Marco Barden, Eric Bell, Ivo Labbe, Daniel 
H. Mcintosh, Hans-Walter Rix and Gregory Rudnick for 
interest and useful discussion. We also would like to thank 
to David Butler for kindly proof-reading parts of this paper. 
We wish to thank the anonymous referee for the helpful 
comments. 

Based on observations made with the NASA/ESA Hub- 
ble Space Telescope, obtained from the data archive at the 
Space Telescope Institute. STScI. STScI is operated by the 
association of Universities for Research in Astronomy, Inc. 
under the NASA contract NAS 5-26555. 

Funding for the creation and distribution of the SDSS 
Archive has been provided by the Alfred P. Sloan Founda- 
tion, the Participating Institutions, the National Aeronau- 
tics and Space Administration, the National Science Foun- 
dation, the US Department of Energy, the Japanese Mon- 
bukagakusho, and the Max-Planck Society. The SDSS Web 
site is [http://w ww.sdss.org The SDSS is managed by the 
Astrophysical Research Consortium (ARC) for the Partic- 
ipating Institutions. The Participating Institutions are the 
University of Chicago, Fermilab, the Institute for Advanced 
Study, the Japan Participation Group, Johns Hopkins Uni- 
versity, Los Alamos National Laboratory, the Max-Planck- 
Institut fur Astronomie (MPIA), the Max-Planck-Institut 
fur Astrophysik (MPA), New Mexico State University, Uni- 
versity of Pittsburgh, Princeton University, the US Naval 
Observatory, and the University of Washington. 



REFERENCES 

Abraham, R. G, van den Bergh S., Glazebrook K., Ellis R. 
S., Santiago B. X., Surma P., Griffiths R. E., 1996, ApJS, 
107, 1 

Abramowitz M., & Stegun I., 1964, Handbook of Mathe- 
matical Functions. Dover, New York. 
Aguerri J.A.L. & Trujillo I., 2002, MNRAS, 333, 633 
Aguerri J.A.L. , Iglesias-Paramo J., Vflchez J. M. & Munoz- 
Tunon C, 2004, AJ, 127, 1344 

Aguerri J.A.L. & Trujillo I., 2004, in preparation (Paper 
II) 

Andredakis Y. C, Peletier R., Balcells M., 1995, MNRAS, 
275, 874 

Blanton M.R., et al., 2003, ApJ, 592, 819 
Barden M., et al., 2004, in preparation 
Bell E.F., Wolf C, Meisenheimer K., Rix H-W., Borch A., 
Dye S., Kleinheinrich M., Mcintosh D.H., 2004, ApJ, sub- 
mitted, astro-ph/0303394 
Binggeli B., & Jerjen, H., 1998, A&A, 333, 17 
Bouwens R.J., Cayon L., & Silk J., 1997, ApJ, 489, L21 
Bouwens R.J., & Silk J., 2002, ApJ, 568, 522 
Caon N., Capaccioli M. & Rampazo R., 1990, A&AS, 86, 
429 

Caon N., Capaccioli M. & D'Onofrio M., 1994, A&AS, 106, 
199 

Cayon L., Silk J., & Chariot S., 1996, 467, L53 

Casertano S. et al. 2000, AJ, 120, 2747 

Chiosi C, & Carraro G, 2002, MNRAS, 335, 335 



16 /. Trujillo and J.A.L. Aguerri 



Coleman G. D., Wu C.-C, Weedman D. W., 1980, ApJS, 
43, 393 

Daddi E. et al. 2003, ApJ, 588, 50 

Eggen O. J., Lynden-Bell D. & Sandage, A. R., 1962, ApJ, 
136, 748 

Fasano G., & Franceschini A., 1987, MNRAS, 225, 155 
Fasano G., Cristiani S., Arnouts S. & Filippi M., 1998, AJ, 
115, 1400 

Fernandez-Soto A, Lanzetta K. M., & Yahil A., 1999, ApJ, 
513, 34 

Franx M. et al. 2003, ApJ, 587, L79 

Freeman K., 1970, ApJ, 160, 811 

Graham A. W., 2001, AJ, 121, 820 

Graham A. W. et al. 2004, in preparation 

Gutierrez CM., Trujillo I., Aguerri J.A.L., Graham A.W., 

Caon N., 2004, ApJ, 602, 664 

Im M., Griffiths, R.E., Ratnatunga, K.U., & Sarajedini 

V.L. 1996, ApJ, 461, L79 

Im M. et al., 2002, ApJ, 571, 136 

Kauffmann G., Colberg J. M., Diaferio, A, & White, S. D. 
M.,1999, MNRAS, 307, 529 

Kauffmann G., & Chariot, S., 1998, MNRAS, 297, 23 
Labbe I. et al. 2003, AJ, 125, 1107 
Larson R., 1975, MNRAS, 173, 671 
Lilly S., et al. 1998, ApJ, 500, 75 

Mao S., Mo H. J., & White S. D. M, 1998, MNRAS, 297, 
L71 

Marleau F.R., & Simard L., 1998, ApJ, 507, 585 (MS98) 
Mo H. J., Mao S. & White S. D. M., 1998, MNRAS, 295, 
319 

Moth P. & Elston R.J., 2002, AJ, 124, 1886 

Peng C.Y., Ho L.C., Impey CD. & Rix H.W., 2002, AJ, 

124, 266 

Press W. H., Teukolsky S. A., Vetterling W. T. & Flan- 
nery, B. P. 1992, Numerical Recipes (Cambridge: Cam- 
bridge Univ. Press) 

Ravindranath, S. et al. 2004, ApJ, 604, L9 

Roche N., Ratnatunga K., Griffiths R. E., Im M., Nairn A., 

1998, MNRAS, 293, 157 

Rudnick G. et al. 2001, AJ, 122, 2205 

Schade D., Lilly S.J., Crampton D., Hammer F., Le Fevre 
O., & Tresse L., 1995, ApJ, 451, LI 

Schade D., Lilly S.J., Le Fevre O., Hammer F., Crampton 

D., 1996, ApJ, 464, 79 

Schade D. et al. 1999, ApJ, 525, 31 

Sersic J., 1968, Atlas de Galaxias Australes Cordoba: Ob- 
servatorio Astronomico 

Shen S., Mo H. J., White S. D. M., Blanton M. R., Kauff- 
mann C, Voges W., Brinkmann J., Csabai I., 2003, MN- 
RAS, 343, 978 

Simard L. 1998, in ASP Conf. Ser. 145, Astronomical Data 

Analysis Software Systems VII, ed. R. Albrecht, R. N. 

Hook, & H. A. Bushouse (San Francisco: ASP), 108 

Simard L. et al. 1999, ApJ, 519, 563 

Simien F., & de Vaucouleurs, C, 1986, ApJ, 302, 564 

Somerville R. S. & Primack, J. R., 1999, MNRAS, 310, 

1087 

Toomre A., & Toomre J., 1972, ApJ, 178, 623 

Trujillo I., Aguerri, J.A.L. , Gutierrez CM., & Cepa J., 

2001a, AJ, 122, 38 

Trujillo I., Aguerri J.A.L., Cepa J. & Gutierrez, CM., 
2001b, MNRAS, 328, 977 



Trujillo I., Aguerri J.A.L., Cepa J. & Gutierrez, CM., 
2001c, MNRAS, 321, 269 

Trujillo I., Graham A.W., Caon N., 2001, MNRAS, 326, 
869 

Trujillo I., et al. 2004, ApJ, 604, 521 

van den Bergh S., Abraham R. C, Ellis, R. S., Tanvir N. 
R., Santiago, B. X. & Glazebrook, K. C, 1996, AJ, 112, 
359 (vdB96) 

van Dokkum P. G. & Ellis R. S., 2003, ApJ, 592, L53 

van Dokkum PC et al. 2003, ApJ, 587, L83 

Vazdekis A., Trujillo I., & Yamada, Y, 2004, ApJ, 601, 33L 

White S. D. M. & Frenk, C. S., 1991, ApJ, 379, 52 

Williams R. E. et al, 1996, ApJ, 112, 1335 

York D.et al., 2000, AJ, 120, 1579 



